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O ■ I. INTRODUCTION 

o : 

' Recently Belle collaboration has studied the production of pseudoscalar and vector charmonium mesons in electron- 
positron annihilation at y/s = 10.6 GeV p|. The lower bound of the cross section measured at Belle 



=3 
< 



a{e^e ^ ipijc) > 33fb 



' is about an order of magnitude higher, than the theoretical predictions 2\. Some efforts were made to explain this 
. discrepancy. For example, in it was assumed, that some of the Belle's 'ipVc signal could actually be iptp events and 

the value of the cross-section of two photon process e+e~ — > was presented. However in later works 0,0 it was 
, shown that QCD corrections decrease this cross section and the subsequent experimental analysis excluded this 

possibility completely. Another possible exclamations could be the contributions of glueball [3 and color octet states 
[ to this process 

' Lately the discrepancy between theoretical predictions and experimental value for the cross section of the process 
I e'^e^ — > iprjc has found the surprisingly simple explanation. In the works 0, it was shown, that taking into 
consideration intrinsic motion of quarks inside mesons one can significantly decrease the virtuality of intermediate 
particles and hence raise the value of the cross section. The method of light cone expansion, used in these works, 
is actually an expansion in a small parameter r ^ M/ y/s, where M is the meson mass. It should be mentioned 
however that this method can lead to large uncertainties, caused by poor knowledge of quark distribution functions. 
Oh, As it was shown in ^] one can significantly change the result by changing the expressions of these functions. Another 
^ ^' method, used prior to these works was the expansion in internal velocity v of quarks in a meson. In most of the works 
' considering the process e~^e~ — > i/^yc only the leading term of this expansion was used and the decrease of virtuality of 
internal particles was neglected. The purpose of our paper is to consider of this effect in the framework of amplitude 
^ ■ expansion over the internal velocity v and comparison of our results with light cone approximations. 
• ^ , Besides the reaction e^e~ tprjc there are also other processes in which the intermediate particles propagate with 
' large virtuality when one neglects the internal motion of quarks inside mesons. Among these processes we would like 
^ , to mention the decays of scalar and tensor mesons into two vector ones, i.e. xo,2 — > VV . In this paper we will consider 
- - ' some of these decays, namely Xc0,2 PP, 4'4> and Xw.2 V'V'- The main reasons for this interest are: 

The only known experimental value for such processes [Tol | 

Br(xcO^'/'0) = (1.0 ±0.6) -10-3 (1) 

is much greater than the theoretical predictions made in the assumption that the intrinsic motion of quarks can 
be neglected (in what follows this assumption will be called " (5-approximation" ) [Til IT^ and one could expect 
that taking into account internal motion would increase them. 

As it was mentioned in ^j, the decays of Xbo,2 can be used to study these mesons in high energy experiments 
(for example in the reaction pp XbjX iptpX). The branching ratio of the decay Xbo ~^ V'V' obtained in this 
work 

Br(xbo V'V') = 3-10-^ 
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is however too small to use this mode in studying the properties of Xbo meson. When estimating this branching 
fraction the authors have neglected the internal motion of c quarks in the hard part of the amplitude. As it 
will be shown bellow, taking it into account one can increase the value of this branching fraction by an order of 
magnitude. In this conditions the observation of Xb mesons in ipip mode is a feasible task. 



II. e+e iprjc 

The diagrams contributing to the process e+e~ — > ijjrjc at the leading order in the strong coupling constant as are 
shown on fig^and the corresponding amplitude can be written in the form 0]: 

2 

M - -^v(fc2)7^u(fci)(V' + ?7c|Jp(0)|0). 
If one neglects the internal motion of quarks in mesons, the matrix element of the electromagnetic current is equal to 

where \I'^,,;^(0) are the values of tp and r]c wave functions in the origin. In this approximation the value of the cross 
section of this process is = 1.9 fb, about an order of magnitude smaller, than the experimental result 0. As it 
was mentioned above, a possible reason for this discrepancy is that internal motion of quarks was neglected in the 
hard scattering part of this process. This effect was studied in the framework of light cone expansion |^ |^, i.e. an 
expansion in the parameter r ~ M/y/s ^ 1/3. The calculations show, that one increases greatly the value of the cross 
section and allows to achieve the same values. It should be mentioned, that the method used in these works can lead 
to uncertainties caused by the poor knowledge of the distribution functions. One can obtain almost arbitrary result 
by changing these functions Q . 

Another method used in studying the process e'^e~ iprjc is the expansion of the amplitude in the relative velocity 
V of quarks in a meson. The potential models [l^ give the value = 0.23, so the expansion parameter v ~ 0.5, 
that is certainly larger, than the expansion parameter of the light cone expansion formalism. In our article we intend 
to consider the decrease of virtuality of intermediate quarks and gluons in the framework of the following simple 
model. The contribution of the diagram fig^ to the amplitude of the e+e~ annihilation into two cc pairs and the 
hadronization of these pairs into final mesons ip and rjc is proportional to 

^ _ u{p'/2 + k')rv{p/2-k) 



{{p + p')/2 + k' - kf{{p + p'/2 + k'Y - ml) 



(f 



+ k)-i^,[k-p' - ^+mc\'^^vi^-k' ] , (3) 



where k and k' are the relative momenta of quarks in tp and rjc mesons, u and v are quark and antiquark bispinors 
and l^^ is the leptonic current. If the motion of quark-antiquark pair in the meson is considered as non-relativistic, 
one can expand the expression (jS)) over the relative momenta k and k' . The expression for example is derived 
from the first term of this expansion. Taking the higher terms of this expansion we however encounter the difficulty 
caused by the strong dependence of propagators in eq. © on fc and k' (actually, at large values of this variables the 
series will not converge at all). To avoid this difficulty we will use the exact expressions for the denominators in eq.® 
and neglect internal motion of quarks in the numerator of this expression. In the c.m.s. of the iprjc pair the following 
formulae are valid: 

p = 7i\/{l,0i,/3}, 
p' = 7Af {1,0^,-/3}, 
k = {79||/3,q-L,79||} , 

k' = {-79||/3,qY,7'?||}, 

where M is the mass of the charmonium meson (for simplicity we neglect the difference between vector and pseu- 
doscalar meson masses), /? is the velocity of one of the charmonia (directed along the Oz axes), q and q' are the space 
parts of the relative momenta k and k' in the rest frames of the corresponding mesons, qn and q'^^ are z-projections of 
these vectors and q^ and q^ are their two-dimensional components perpendicular to this axes. From these equations 
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the following expressions for the scalar products can be obtained: 

pk' = -2/372 Afg[|, 
p'k = 2/37'Mg||, 
kk' = -7'9||«||(l + /3')-qiq'±, 

= -gj|-qi, 
= - q± • 

It is clear that the scalar products (pfc') and (p'fc) are enhanced by the factor 7 — ^/s/2M ^ in {kk') there are both 
enhanced and unenhanced terms and scalar products k'^ and k'^ do not contain this factor at all. In the light cone 
expansion formalism the perpendicular components of vectors q and q' (i.e. terms that are not enhanced by 7) 
are neglected. This approximation is valid for reactions with high cm. energy, where 7 3> 1, but can lead to 
sufficient errors in other cases. That is why we will carry out our calculations using two methods — both neglecting 
contributions of this terms and leaving them nonzero and afterwards we will compare the results. 

The internal motion of quarks inside ip meson will be described by setting its quark and antiquark momenta to 
P/2 + fc and P/2 — k respectively. Similarly, the momenta of quark and antiquark hadronizing into r/c will be P'±fc'/2. 
After this substitution the gluon propagator will have the form 



k'-k] = l[l + Z{p'k^pk') + -{k~k'f]=-± (4) 



2 / 4 V s ' ' s / 4 

and quark propagators on fig^will be equal to 

l+p'-k^-ml = i(l-V + ^(fc^ + ^-rn^))=^.„ (5) 



From equations lO, ^ it is clear, that in (5-approximation ( i.e. neglecting internal motion and assuming the 
quark mass = M/2) the dimensionless propagators d and si 2 will be equal to 1. Keeping this in mind it is easy 
to understand that to obtain the expression for any diagram of fig ^ one should divide the expression where 
this motion was neglected, by this dimensionless propagators and integrate it with the proper chosen wave functions. 
It is useful to notice, that the diagrams can be obtained by charge conjugation of the diagram QJl, so it is not 
necessary to calculate it separately (the same is also valid for diagramsCJ) andOJi) Thus, the amplitude of the process 
e^e~ — > Tprjc can be written in the form 

M = A^otV / dp^, f — + — ) = MqK, (7) 



J dp J 2d\si S2 

where the integration is hold over the wave functions of "ip and 77c mesons in momentum representation: 

dp = d^qd^q'(j)^{q)(l)^^{c{). 

In equation Q A^o is the amplitude of this process with no internal motion and the factor K describes the decrease 
of momenta of virtual quarks and gluon. The variation of wave function width results in change of this factor and 
therefore the variation of annihilation cross section. It is clear, that if we fix the value of the amplitude and tend the 
width of the wave function to zero (i.e. neglect the internal motion of quarks), than the factor K will tend to 1 and 
the equation iQ will reduce to (|2Jl. It should be mentioned also, the values ^"7/^(0), ^'r;c(0) ^Iso depend on the choice 
of the wave functions — it is obvious, that with increase of the width of the distribution in q the values of the wave 
function in the origin will grow. Thus, the variation of wave function width will lead to double effect: the variation 
of the momenta of internal particles and the change of the values of wave functions in the origin. Both these effects 
are taken into account in the equation |7I). 



III. WAVE FUNCTIONS 



In our model we assign the momentum P/2 + fc to c quark in ip. It is evident, that this quark is not on mass 
shell: its virtuality depends on the relative momentum k and equals to + k"^ . The formal use of this formula 
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can lead to unphysical situations, when the virtuahty of c quark is zero or negative. It is obvious, that such cases 
would occur only if the motion of cc pair in the meson is relativistic, whereas in the model used here is valid only for 
non-relativistic movement, when the virtuality of this quark does not differ significantly from A//2. The calculations 
described below will show, that even in the broad wave-function case ( = 0.4) the maximum of the amplitude is 
accumulated at — fc^ ~ 1 GeV^, that corresponds to virtuality ^1.2 GeV. 

Another unphysical singularity caused by the fact that the virtuality of quarks and antiquarks can be zero or 
negative is that the propagators l@J, (O,© can vanish at some values of quarks' momenta. As a result one can 
obtain almost arbitrary values for the cross section of the process under consideration by using different distribution 
functions. 

In the previous works the potential quark model was often used to derive the form of the wave function. Such a 
function, obtained from Schrodinger equation solution allows, although with a small probability, any values of the 
relative momentum of quarks and we will get an infinitive result when substituting it in the equation ( 01 • This 
problem is obviously caused by the singularity of propagators in the region q ~ M /2, where the motion of the quarks 
is relativistic and the usage of non-relativistic wave functions is not valid. To avoid this difficulty, when deriving the 
expression for the wave function we will use another method, proposed in Ti]. 

In what follows we will restrict ourself to the valence quark approximation, i.e. we will suppose that t/i and rjc 
mesons are built only from cc pair. Let the momentum of the quark in the meson rest frame he q — (goiQ) E^nd the 
meson's wave function in the momentum representation is ^(q). Then, in the frame of infinitive momentum the c 
quark distribution function will be 



M 

where M is the meson mass and x is the c-quark momentum fraction. For the mesons with hidden flavour in valence 
quark approximation we have go = A//2. It is useful to introduce the following integration variables: 

a q ~ — - — vdvdq\\^ 



where v — 2|q|/M. After integration over q^ 



/(-) - ^ /' dv^\<i>{v^)\\ 
= (l-2a;)^ 

Differentiating this relation with respect to x we get 

From the last formula one can obtain the relation between meson wave function and the structure function in the 
infinitive momentum frame: 

Let us use the Regge parametrization of the structure function 0, 0, 0| : 

Since we restrict ourself to two-body state cc and neglect the contribution of the sea quarks the parameter 7 should 
be set equal to zero. Substituting the expression for the wave function into (jS]) we get 

\ 1/2 /I 2\ (-c«-l)/2 



Let us determine the value of the parameter a. We can do this by fixing either the parameters of the V' r^gge 
trajectory, or mean square of the velocity of c-quark in meson. In what follows we will use the latter method. The 
calculations will be held using four different values of mean quark velocity :u^ = 0.2, 0.3, 1/3 and 0.4. These values 
correspond to a = —6.0, —3.5, —3.0 and —2.25. It is interesting to notice, that for — 1/3, a = —3.0 0| the wave 
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function integrated over transvese momentum differs from light cone function used in Q no more than by 5%. We 
have chosen the value — 1/3 because the calculation of ^ e^e^ and ?7c ~* 77 decay widths with this function 
leads to the values r('0 e+e^) = 5.45 keV, r(77c — > 77) = 7.3 keV, that are in agreement with the experiment 

Elil. 

The results of our calculations for listed values of are presented in table HJ In the first column there are values 
for |i?(0)p = 47r|^'(0)p, the second one contains the values of the cross section obtained from equation ^ neglecting 
internal quark motion in the hard part of the amplitude, the last two ones contain the values of the factor from 
formula Q neglecting and leaving the transverse components of internal momenta k and k' . Finally, in last column 
we give the values of the total cross section. 

The results of our calculations, presented in table show, that the calculation held with exact dependence of 
propagators on internal momenta gives the values that are lager, than the values obtained neglecting this motion. 
For example, for the case = 0.4 the cross section is multiplied by 4. So our assumption that the dependence of 
propagators on the internal momenta is strong proved to be correct and the expansion of this propagators in relative 
velocity v leads to the series with rather poor convergence. Although, the effect described above increases the value of 
the cross section, this enhancement is not large enough to explains the experimental data obtained at Belle. As it was 
mentioned above, the calculations in the framework of the light cone expansion method ^] give the values that are in 
agreement with experimental data and, if there is no other reasons that can increase the cross section obtained in the 
framework of NRQCD by an order of magnitude we can assume that the light cone expansion parameter Mj ^/s ~ 1/3 
is more suitable for describing the double quarkonia states in the electron- positron annihilation at — 10.6 GeV, 
then the relative velocity v ^ 0.5. 

Another interesting feature that can be seen from table ^ is that taking into account the terms suppressed by s 
leads to significant variation of the results. For the narrow wave function with = 0.2 this correction is about 
20%, while for the wide one with v'^ = 0.4 the cross section is almost doubled. Thus, the transverse components of 
internal momenta have to be taken into account. In the light cone expansion formalism the terms suppressed by s 
are neglected and it results in a large error. However, taking this terms into account will only increase the one, and 
a new result does not contradict the experiment. 



IV. XJ 



vv 



In addition to the 



^rjc process considered in the previous sections there are also other reactions in which 



intermediate particles have large virtuality when (^-approximation is used. Among these processes we would like to 
note the decays xo,2 VV. In the (5-approximation these decays were already studied in literature [Til ll^ 

Two diagrams contributing to such decays at the leading order on the strong coupling constant ag are shown in 
fig[21 the others can be obtained by interchanging final vector mesons. It can be easily seen, that the virtuality of 
gluons on these diagrams in the (5-approximation equals M^/4 and one could expect, that similarly to the results of 
previous sections it will decrease if the internal motion of quarks in the mesons is taken into account. 

Analytical formulae for the width of the decay xo ~* VV were presented in the work '20^ and we will use this results 
in our paper. The nonzero helicity amplitudes A^^^ of the decay of scalar meson xo into vector mesons Vi and V2 
with the helicities Ai and A2 are given by the expressions 



•^1,1 



-^0,0 



A 



(0) 



)13 



-i — ^TT^a? 



9V3 
i?'(0)| 



\R'm 

M4 



M 



where 



m/M, 



m and AI are masses of vector and scalar mesons respectively, R{r) is the radial part of the scalar meson's wave 

function, /l,t are longitudinal and transverse leptonic constants of the vector meson and coefficients are equal 

to ' 



AO) _ 
-'1,1 — 



1 

-— / dxdy(i3T{x)(l)T{y)- 
0/ J 



1 



1 



xy + {x — y)^£^ (1 — a;)(l ~ y) + (x — y)^e^ 



2xy — X ~ y + 2(x — y)^e^ 



1 

1 + :^ 



2 2a;?; — x ~ y + 2{x — y)^e^ 
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7-(0) 

^0,0 







1 



1 



1 



(x-j;)2(l-4e2) 



2xy — X — y + 2{x — y)^e^ \ 2 2xy — x — y + 2{x — y)^e^ 
1 (x-y)2(l-4e2) 



-2e' 



1 



2 2xy — a; — y + 2{x — j/)^e^ 



In the above equations x and y are the momentum fractions of the final mesons, carried by quarks and 4>l,t{x) are 
longitudinal and transverse distribution functions of these quarks in mesons. Later we will show, that the contribution 
of longitudinally polarized mesons to this process is small and one can safely neglect it, as it was done in [TH . 
In | 20j| the similar formulae for nonzero helicity amplitudes of tensor meson decay are also presented: 



•^AiAaiAi ~ -^AiAa^ "m,Ai -Ai l^^J ' 

where /z is the meson spin projection on fixed axe, 9 and are polar and azimuthal angles of one of the final mesons 
in X2 rest frame and reduced amplitudes Ax^x^ are given by the expressions 



{2) 





^1,1 








= A,^,=A-,^,^A,^-,= ^''/'-^a^.'^^'^'/T/./l^^, 












.2"y2 |i?'(0)| 2^(2) 


r(2) 
1,1 


_ r(0) 






1 

^32 


1 1 

/ / dxdy(j>T{x)(f>L{y) , ,2 2 n ^n ^ , ^ \2 2 
J J xy + [x - yye^ [1 - x)(l - y) + [x - yye-^ 




7-(2) 
1,0 



2xy — X — y + 2{x — y)^e^ 



1 



(x-y)2(l~4e2) 



2 2a;?; — a: — ?/ + 2(a: — y)^e^ 



^^'^11 dxdy(l)T{x)(l)L{y) 



1 



xy + {x — y)^e^ (1 — — y) + (a; — j/)^e^ 



2x?/ — a; — y + 2(a; — yYe^ ' 



1 1 



4j] = y dxdy(t)Tix)(t>L{y) 





2xy — a: — y + 2(a: — y)^e^ 



1 



1 



a;y + (a; — y)^e^ (1 — a;)(l — y) + (x — y)^e^ 

(^_y)2(l_4,2) 



1 



2xy — X — y + 2{x — y)^e^ 



1 

1 + ::: 



(a; - y)2(l - 4e2) 



2 2a;y — a; — y + 2(a; — y)^e^ 



In the literature plenty of different expressions for distribution functions (t>L.T{x) can be found (review is given in 
|2l|) and, as it was shown in [sj, the result depends strongly on the their choice. That is why in our calculations we 
will use several variants: 
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• Chernyak, Zhitnitsky ^2d\ (in the tables this set will be labeled as "CZ") : 

(j)L (x) = 0(^2 ' (x) = (j>2 (x) , 

• Assymptotic ^ ("0i") : 

(j)L{x) = </)T (x) = (/)! (x) , (10) 

• Symmetric [3 Ch") 

(j)L {x) = (/)T (x) = (/)3 (x) , 

• (^-approximation ("(5") 



!)L (x) = 0T (a;) = (5 I X - ^ 



All these distributions are normalized by the condition 

1 

jHx)dx = 1 



and the effective width of the distribution 6 is defined as 







In the above equations 



0CZ {x) = 13. 2x (1 - x) - 36x2 _ 


x) , 




w 6.7-10-2 


4>i {x) = 6x(l 


- x), 


(^^>,, « 


^ 5 - 10-2, 


02 (x) = 30x2(1 - 


'X)\ 




i 3.6 - 10-2, 


03 (x) = 14Gx^ (1 - 


x) , 




i 2.8- 10-2. 



By analogy with the reaction e^e~ — > ■(/"7c one can expect, the branching fraction of the decay Xbo,2 V'V' will 
increase with the increase of the width of the distribution. 

In our work we consider the decays XcQ.2 PP, XcO.2 — > 00 and X60,2 — > '0V'- Numerical values for parameters are 
listed in table ^ The values of branching fractions and decay widths obtained with the help of these distributions 
are presented in table IIIII From this table one can easily see that the branching fractions do strongly depend on 
the distribution functions and increase with increase of their "widths". As it was mentioned above, the experimental 
value of the Xco — > 00 branching fraction is much greater than the theoretical predictions made in the framework of 
(^-approximation. Now internal motion of quarks in mesons raises these predictions. For example, with the distribution 
JTUJl we get 

Br(xco^00) = 0.9-10-3, 

while the experimental result is (1.0±G.6) - IQ-"^. The same can be said about the decays X6o,2 V'V'- The calculations 
held in (5-approximation give Br(xf,o V'V') = 1-4 - 10"^, and now this value can be raised to 2 - 10-"*. 

In table IIVI we present the ratio of helicity amplitudes for production of longitudinally and transversely polarized 
vector mesons. From this table it is clear that this ratio is less than 15% for all the considered decays and the 
production of longitudinally polarized mesons can be safely neglected, as it was done in [T]| . 
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V. pp ^l)^f}X 

As it was mentioned in the decay Xfco,2 V'V' can be utilized in experiment to study the properties of xbJ 
mesons in high energy reactions (for example, in the inclusive process pp ggX Xbo,2X ~^ ipipX). It should 
be mentioned, that in addition to resonance production of ^''4' ps^ir via Xbj decays there is also a non-resonance 
production of tjjijj pair in continuum and it is important to study the possibility of separating of the resonance signal 
from this background. 

The cross section of the double tp production in the proton-antiproton annihilation can be related to the cross 
section of its partonic subprocess a{ab — > ipip) and the structure functions of the initial hadrons according to the 
relation 

a{pp ^ abX ^ i^i^X) =Y.j d^ad(,,f^\^a)fi^\ib)o{ah ^ U), (H) 

where and are the momentum fractions of partons a and b, fa^\^a) and fl,^\£,b) are the distribution functions 
of partons a and b in proton and antiproton. It is convenient to use new integration variables s w s^aCf) — the square 
of the invariant mass of the partonic pair and ^ = — ^b- If we are interested in the central production of ipip pair 
in the resonance region, we have = ^b = -^x/v^ 1- this case the gluonic pair (i.e. a = b = g) gives the main 
contribution to the reaction give and the equation can be written in the form 



da{pp TpipX) _ 1 o-gg 
ds 16 s 



L, (12) 



where the first factor corresponds to averaging over gluon colors and helicities, agg = a{gg i^i^), and the dimen- 
sionless factor 

l-s/s 

m = 2 j J^j(p)i^,)f(p)i^,)d^ (13) 



describes the partonic luminosity. 

The partonic cross section agg equals the sum of cross sections of resonance reaction pp — > ggX — > XbJ^ ~^ ^pipX 
and the non-resonance background: 



Resonance part of the cross section can be obtained from the values presented in table IIIII with the help of Breit- 
Wigner formula: 

= M^(^'^ + ^) (M ~M )2 , r2 /4 ' 

,7=0,2 93 ^ ^"Xb.ll +^Xbjl* 

where Mgg = ^fi — invariant mass of the colliding gluons and total widths of x&j-mesons can be obtained from the 
values of the radial part R{r) of their wave functions in the origin p^ : 



rx.0 = nXbo - gg) = 96-^|i?'(0)|2 « 0.544MeV, 

XbO 

^Xb. = r(x62^55) = ^T|f-|i?'(0)P«0.14MeV. 

Xbl 
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The implementation error A could however be much smaller than this widths and it can be useful to take it into 
account. This can be done by means of simple substitution 

It can be easily seen that this substitution does not change the value of the total integrated cross section (i.e. 

JardMgg= Ja^dMgg). 
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Non-resonance part of the cross section (i.e. cross section of the double ip production in continuum) was obtained 
in I23, 01 • According to this works 



4 

S4 



^4 ^ y"C', 

dcT 3 
di ~ 47r 



\ i ) \ ) (t^Miy{u-MlY' ^ ^ 

where s, t and ii are the Mandelstam variables of the partonic subprocess, Af^ and V'(O) are mass and wave function 



of V', 2/ = {i + u)/2s 


, z = (^- 


Co = 


^(335- 


Ci = 






-(-1- 
3 ^ 


C3 - 


64 , 


C4 = 


256 (1 + 



^ 13568 + 98144 z"^ + 35584 - 1186048 z^ + 258048 2^° + 3981312 z^^) ^ 
' 14288 z^ + 168064 z* - 403840 z^ - 596224 z^ + 1935360 z^") , 
11472 z^ + 239008 z^ - 1010944 z^ + 1134336 z^) , 
828z2 + 688 z^ - 10176 z^), 

8z^ + mz^) . 

On fig|3|the dependence of partonic cross section form the invariant mass of colliding gluons in the resonance mass 
region and the value for the implementation error A = 50 McV is shown. It is clearly seen from this figure that the 
separation of Xb signal from the background for most of distributions should not make any troubles. On the other 
hand for the (^-approximation the signal/background ratio 



= 



a. 



A 



is not high enough and it may be useful to study the possibilities of increasing this ratio. It can be made by applying 
certain cuts on the kinematics of final ■(/' mesons. As it can bee seen in fig01 where the angular distribution of the 
background cross section over the cosine of scattering angle x = cos 9 is shown, final ip mesons fly mainly along 
the beam axes, while the distribution of "0 produced in Xb decays is isotropic. Thats why applying the condition 
\x\ < Xmax we Can increase the ratio Ti/^ . On fig,|31the dependence of this ratio form the cut angle is shown. 

Let us now return to the total cross section of the Xf,o,2-nieson production. According to formula ifT^ this cross 
section is equal to 

E ^{PP^Xb.J + X) ^ f $L(s)a.(l), (17) 

J=0,2 * 

where ar{s) is the resonance cross section ()14|l . Because the resonances are narrow, we can neglect the dependence 
of the partonic luminosity L{s) on s, so the integral over s in H17|) will reduce to a simple expression: 

<y{pP XbJ + X) = -^L E (2J+l)r(x6j^55), 

J=0,2 X J=0,2 

where the value of the partonic luminosity L is defined according to formula l|13j) . For the gluon distribution functions 
in the initial (anti)protons we have used the parameterization presented in [25j at the value of virtuality — Af^. 
The values of Xf,o and Xb2 production cross sections obtained using this formula are equal to 

ct(pp ^ Xho + X) = 250 nb, a{pp Xb2 + X) = 320 nb 

at the Tevatron energies (^/s = 2 TeV) and 

a{pp -> Xbo + X) = 1.5 ph, a{pp ^ Xb2 + X) = 2 /ib 

at LHC energies {y/s ~ 14 TeV). In the tpip mode these values correspond to 

a{pp -> XbQ + -'^)Br(x6o H^) = 53 pb, a{pp -> Xb2 + ^)Br(xb2 ^ V'V') = 170 pb 
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at = 2 TeV and 

<^{PP Xbo + -'^)Br(x6o i'i^) = 330 pb, a{pp Xb2 + ^)Br(x62 V"/") = 1 nb 
at = 14 TeV. 

The possibility of quarkonia observation in their decays was also proposed in the works . The value of 

the r]b production cross section at Tevatron presented there 

a{pp ijb + X) = 2.5 fih 

is about an order of magnitude higher, then our results. This difference can be explained by the difference in the total 
widths of r]b and Xb mesons. Inspired by these works the CDF collaboraion has measured the ^/j^ mass distribution 
and has observed a small cluster in the r]f, region ^. Since this cluster is small only upper 95% CL upper limit 
on cr,,^(|y| < 0.4)Br(77b 'lp'^p)[Tir{^p /i"''^")]^ was found in this experiment. As it can be seen from the ipip mass 
spectrum, there are also some events in the Xb region, so one can hope that with the increase of luminocity it will 
become possible to observe Xb mesons in the reaction studied in this paper. 

It should be mentioned, that our estimates were obtained in the coUinear gluon approximation, i.e. using the 
integrated over the transverse momentum gluon structure functions. As a result of this approximation we have 
underestimated the values of the Xb meson production cross section and it is comparable with the cross section of 
direct T production [2^. So, our estimates should be considered as lower bounds and the situation could be more 
favorable for the Xb observation in tp^jj mode. 



VI. CONCLUSION 



The discrepancy of experimental results and theoretical predictions for the cross section of production in e^e~ 
annihilation have found a natural explamation by taking into account intrinsic motion of quarks in ip and ijc mesons 
in the amplitude of heavy quarks production. At this point our results coincide with the results of the works . 

We have also studied another example of the reaction xo,2 VV, and specifically the decays of Xfco(0''"'') and 
Xb2(2^^) into a pair of ijj mesons. It is shown, that taking into account the relative motion of quarks in the amplitude 
of the decay of Xb meson to c-quarks increase the branching fractions of these decays by an order of magnitude. The 
agreement with the experiment was achieved for the only experimentally known decay XcO ~^ <t'4>- 

Later we have considered the possibility of the observation of Xb V'V' decay in the experiment. The values of the 
cross sections of Xb production in pp annihilation and its subsequent decay into iptp pair presented in our paper show, 
that it is possible to observe this reaction at Tevatron and LHC energies. 
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TABLE L Cross sections of the e^e — > i/"?c process using different parametrizations of meson wave functions 



«2 


|i?(0)p(GeV) 


<To(fb) 


cr/o-f) 




a(fb) 


0.2 


0.31 


0.78 


1.9 


2.3 


1.7 


0.3 


0.50 


1.9 


2.4 


3.5 


6.5 


1/3 


0.55 


2.3 


2.6 


4.0 


9.2 


0.4 


0.65 


3.2 


2.9 


5.2 


16.6 




FIG. 2: XJ ^^^^ 



TABLE II: Parameters 





M, 


m,GeV 


Cts 


/T,GeV 


/i,GeV 


|7?'(0)p,GeV^ 




3.415 


0.776 


0.3 


0.2 


0.2 


0.11 


XcO ^ 


3.415 


1.02 


0.3 


0.2 


0.2 


0.11 


X60 ^V' 


9.8599 


3.09692 


0.2 


0.4 


0.4 


1.34 


Xc2 PP 


3.55626 


0.776 


0.3 


0.2 


0.2 


0.11 


Xc2 ^ 0<?!> 


3.55626 


1.02 


0.3 


0.2 


0.2 


0.11 


XM ^ V'V' 


9.9126 


3.09692 


0.2 


0.4 


0.4 


1.34 



TABLE III: r(x ^ VV), keV / Br(x -> VV), 10"'' 





cz 




03 


5 


XcO 


PP 


24.9/35.6 


8.98/12.9 


2.75/3.93 


1.36/1.94 


XcO 




16.3/23.3 


6.31/9.02 


2.10/3.01 


1.09/1.57 


Xbo 




0.117/2.15 


0.0462/0.849 


0.0157/0.288 


0.00827/0.152 


Xc2 


PP 


6.94/43.8 


10.6/66.6 


3.39/21.4 


1.89/12. 


Xc2 




7.14/45.1 


9.79/61.8 


3.38/21.3 


1.96/12.4 


X62 




0.0757/5.32 


0.100/7.05 


0.0355/2.5 


0.0208/1.47 



TABLE IV: 2|yt'^,l/X"^i^% 





CZ 


4>i 


03 


S 


XcO PP 


0.306 


1.99 


2.18 


2.65 


XcO -> H 


1.23 


7.41 


7.95 


9.43 


Xbo ipip 


1.61 


9.59 


10.3 


12.1 



a, pb 

3r 




FIG. 3: Distribution of the partonic cross section versus gluon-gluon mass 




FIG. 5: Signal/background ratio for different distribution functions 



